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ABSTRACT 


A digital simulation language for the interactive definition and 
solution of piece-wise continuous non-linear ordinary piigeeenial 
equations using the Runge-Kutta method has been designed and imple- 
mented. The combination of interactive eae approach and a 
special differential equation description language make the analysis 
program very versatile and easy to use. For second order systems, 
a grid of phase trajectory segments over the user specified phase 


plane is used as a background for the solutions. 








II. 


Uae 


TABLE OF CONTENTS 


onion. 2 09 .. 2... Ses. 8 
Chin RAPE SeGRIrPPrioN OF THE PROGRA Messse2-_----- 13 
A. INTRODUCTION ----------------------------------- 13 
Bo hl eee TiCALD BASIS OF THE PROGRAM ------.--- 14 
Og CT is Ca 16 
D. Bbeeret ONPR HE COMPILER --..-----.ss2---<2 18 
Ee, DES@Rer TIO hrOReGHE INTBRPRE TOR .....-.-----.- 20 
EF. DESCRIPTION OR THE SLOPE LINE*GENERATOR-..- 23 
G. DESCRIPTION OF THE SOLUTION GENERATOR ------ 24 
GENERAL DESCRIPTION OF THE LANGUAGE ------------ 25 
A, INTRODUCTION TO THE LANGUAGE ----------------- 25 
Bes eee oO ht Pipe eNGUAGE ----==-- ----------=-=- 25 
|, Mem@liawaetier Set - - si - cee - - ie sie - - ee - = 25 
ne. = = ae Se ee ae a = eee 26 

De eee Goose a eS So eee oe oe oe ie 26 

4, Expressions ----------------------------------- a 

er Seecliaep C MH xpressions —---=----=-----— = — — — 2a 

be loepical Hxpressions --.---------.-.--=.-_-... 28 
Sweriatwetnatical Functions —-<-.<--<..---<2ss--<s---..- 30 

6. feiwametic Assignment Statements -----..-.-..--.-- 31 

TO @arrOl Statements .«2.2l2.ee26c256cco5eaeoe5 22545 Se5= 33 
Ceeene Shoe Statement <-—..2e.-. = see eee eek 33 

pee tne Arithmetic IF Statement 2-4. 2-5-2... 33 

ee Wieterential Equation Statements -2....=-.<-------<- 33 





G. USING THE LANGUAGE --------------------.-....-- 35 


1. General Comments ---------------------------- 25 
2a a oesiele Witterential Equation ------------------ 35 
3. More than One Differential Equation ------------- 36 
4. Using the Remaining Features of the Language---- 37 
ee en Se METHODS OF USE -------=----=------ a 
A. GENERAL COMMENTS ----------------------------- 39 
B. GENERAL PRINCIPLES ---------------------------- 39 
Pepe yetemm Description -------<-<---------.-----<.- 39 
Zyeiscesewminte Constants -----------~--s------- +--+ 39 
Zee hese Plane Window =-=------2.2---.-.-.-.--- 40 
A ttOnmiGtatMete rs <-.=|=---<-.--2-<-555255se5e 4] 
i onoeltG DETAILS OF THE BATCH VERSION .-----. 41 
Mee ioroetOrapwe Batch Run ------<5---.-sss-<22 41 
Zee @omduAnanancement ----.-.--..-.---.<----- 42 
boa elem bata Deck -..--si68dssccanseoccecacsee 43 

D. SPECIFIC DETAILS OF THE INTERACTIVE 
Cre weteeo VMiinvohON .--.----..--<- je oe es ee eo secee a4 
l. Planning for the Interactive Session’ --..-.-.-.-.- 44 
2. Check-off List for Loading the Program --------- 44 
3. Interacting with the Program ------------------- 46 
4. Suggested First Session Example --------------- 47 
Vv. APPLICATIONS TO VARIOUS PROBLEMS -------------- >t 
ere ere OnmininiIh>S <---------==-=======—-=—--=- 51 
B. VARIOUS SINGLE EQUATION EXAMPLES ----------- 52 
C. VARIOUS TYPES OF PENDULUM EXAMPLES ------- 66 
Pao wee PE) oom bis ~--------------------------- (4 
EIN MSS IG) -----—--—- --— == ——- = a 76 


ad oa | 
= 
= eee 





ep er ee -—----------------------..-.---.-___--- 86 
Pe WW eeND ORDER LINEAR SYSTEM ----------------- 90 
Pee Wu hi @emER SYSTEM -------------------------- 95 
i) HvameRHsisc AND DEAD ZONE ----------- Perea 97 
fe i io Ure S TORING FORGES ------------------ 99 
eer ete Seis VST Mio -----------------+-225------- 106 


VI. CONCLUSIONS AND SUGGESTIONS FOR FURTHER WORK 


IN THIS AREA ----------------------------------------- 110 
APPENDIX A: Required Job Control Language for Batch 
Version ----------- 9-9 oer ee ne rrr eee We 
APPENDIX B: Required Control Cards for Interactive 
Graphics Version ---------rssse- cern eer ec ceene Lis 
PROGRAMA Listing of Batch Version --------------2-------- 114 
PROGRAM B Listing of Interactive Graphics Version ---------- 154 
NS a) ee ieee ee oc eeeiee- --- --2 set 218 
pee) See © NebIS T aman ~~ > oot 219 
FORM DD 1473 --------- rr rere rr rr rrr rr rrr cer rccerere 220 





Figure 
2-1 


2-2 


5-4 
5-5 


5-6 


5-8 


5-9 


Page 
Overall block diagram -------------------------------- ie, 
Compiler simplified block diagram --------------------- 19 
Interpreter ------------------------------------------ 21 
Example 1: X!' - A*(1-X**2)*X' + B¥X = 0 ------------- 53 
Example 2: X‘'' + A¥X*X' + B¥X = 0 eel 55 
Example 3: X!'' + (1.0 - X**2)*X' + X = 0 --------------- 57 
Example 4: X!'' + X'**#2 + X= 0 -e------------2--------- ao 
Example 5: X'' + A*X + B¥X**3 = 0 ------------------- 61 
Example 6: X'' + (1.0 - ABS(X))*X' + X = 0 ------------ 63 
Becamupre (/; X'' + X' + X*ABS(X) = 0 ------------------ 65 
Example 8: X''+X=0 ------------------------------ 67 
Example 9: X'!' + SIN(X) = 0 -------------------------- 69 
Example 10: X'' + A*¥X + SIN(X) = 0 ------------------- 71 
Example 11: X'' + A*X'*ABS(X') + SIN(X) = 0 ---------- (ee 
Example 12: Saturated servo system ---<--------------- 75 
Example 13: Ideal relay ----------------------------- 77 
Hewmple 14: X''+ 0. 2*X' +1. O*SIGN(4. OFX-X')=0. 0 ----- tg? 
Peale 15: X''+ 0.2*X'+ 1. OFSIEGN(4. 04X47 xX')=-0.0 ----- 81 
emiple 16: X''+ 0,2*X' +1. O*SIGN(Z. OF XTX')=0. 0 ----- 83 
Example 17: X''+ 0.2*X'+1. O*SIGN(2. 0*X-X')=0.0 ----- 85 
Example 18: Dead zone ----------------------r-r- cre 87 
Example 19: Dead zone -----------<2--------- 2-5 ----- 89 


LIST OF DRAWINGS 


Example 20(a): Second order linear system,damping=1. 0--94 


Example 


20(b): Second order linear system, damping=0. 5--93 





5 -23(a) 


5 -23(b) 


5-24(a) 


5-24(b) 


5 -24(c) 


Example 
Example 
Example 


Example 


Example 


Example 


Example 
Example 
Example 


Example 


20(c): Second order linear system, damping=0. 1---- 94 
21: A fifth order system ------------------------ 96 
22: Relay with dead zone and hysteresis --------- 98 
23(a): X'' + 0. 9%X! + (KX + AXK#53) = 0.0 

A = -0.1 --e--- renee ene t------------- 100 
23(b): A= 0,2 ------------------------------- 101 
24(a): X'' + 0.9*X' + (X + 0. 2*X*%3 )=A*COS(BFT ) 

Peewee, B90. 7 ise ee 103 
Oe IN SOS) eo ae ee eee 104 
24(c): A=5.0, B=1.0 ------------------------ NGS 
25: A second order bang-bang system------------ 107 
26: A third order bang-bang system ------------- 109 





I INTRODUCTION 


The goal of this study was to develop an interactive analysis tool 
for the study of piecewise continuous non-linear ordinary differential 
equations. It is recognized that there are available in most computer 
centers several programs for solving non-linear differential equations. 
However they are either very special purpose or else require that the 
problem and/or control description be inserted wie a subroutine or 
partof the analysis program. This in turn must be compiled by the 
computer before execution time. If the user wants to change some 
feature of the system, he must change the applicable program cards 
and recompile. Thus in general only one problem type per computer 
run is possible. 

Also there exists a large number of continuous digital simulation 
languages [Clancy and Fineberg, 1965]. These are best represented 
at this time [Cardenas and Karplus, 1970] by MIMIC, DSL/90, and 
CSMP. DSL/90 and CSMP, which is based on DSL/90, have transla- 
tors which convert the prescribed input language into FORTRAN sub- 
routines. These subroutines are then compiled by the FORTRAN 
compiler and executed with the solution routine. This is a multi-step 
process which is both time consuming and not compatible with the 
interactive approach. CSMP and DSL/90 will accept the problem 
description either in the block diagram form or the differential equa- 
tion form. There is a version of CSMP (similar to the above CSMP 
in name only) which uses the block diagram approach only and is 


suitable for use on the small computers of the IBM line. MIMIC 


incorporates its own compiler and has a very flexible algebraic 
capability. 

The program described in this thesis uses the differential equation 
approach and has a built-in compiler flexibility similar to MIMIC. In 
addition it has been designed from its conception to be fully interactive 
using graphics as the interface between the user and the program. All 
problem description is handled as data. Thus only the compiled version 
of the analysis program is needed by the user. This is accomplished 
by the use of a special differential equation description language. An 
automatic compiler (transparent to the user) built into the program 
provides the interface between the description language and the fourth 
order Runge-Kutta solver. This use ofa special differential equation 
description language and the interactive graphics approach give the 
program considerable flexibility in the type of differential equations 
it can handle. 

The class of piece-wise continuous non-linear ordinary differential 
equations handled by this program is restricted by the following con- 
siderations: 

1. The equation must have only one independent variable. 

2. The equation must be normalized so that the coefficient and 
exponent of the highest order derivative are both one. 

3. At each given moment during the solution one and only one 
differential equation can be applicable. The applicable diffe rential 
equation must be well defined for the conditions at that point. Note 
that this allows piece-wise continuous differential equations with 


discontinuities but excludes simultaneous differential equations. 


7. 





These are the only restrictions inherent in the method of solution. 
The implementation described in this thesis has been arbitrarily 
limited to differential equations of eighth order or less. Fora given 
problem there may be added restrictions in defining the differential 
equations due to the limitations of the program language used to define 
the systems and/or the ingenuity of the user. However, if the user 
can define his system using the defined programming language and the 
above three requirements are not violated then the program will solve 
the system. 

Although the language used in the system description is quite 
natural to most users, it is important to completely specify it. Thus 
Chapter III deals with a description of the elements of the language 
and the principles of its use. 

Other features of the program are: 

1. The program has the generality required to be useful for a 
wide variety of applications. In addition, the program does not re- 
quire sophisticated computer techniques on the part of the user. 

2. The outputs of the program are a smooth x vs x' plot (phase 
plane in the second order case) of the solution and printer plots of 
the time response of each state variable. 

3. The ''window'"' represented by the x vs x' plane can be adjusted 
in size, scale, and position. These three parameters can be adjusted 
in the x and x' directions independently; thus the user has full flexibil- 
ity to adjust the window as desired. 

4. The solution is accomplished using the state variable approach 
with fourth order Runge-Kutta integration. The parameters of the 


solution are completely under the control of the user. 
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5. For second order time independent systems there is super- 
imposed on the phase plane a grid of phase trajectory segments. 
These slope markers represent the slope of the trajectory if it were 
to pass through that area of the phase plane. Thus without obtaining 
a solution one can visually obtain a feel for what the system will do 
for various initial conditions. This in turn aids the user in deciding 
what the essential characteristics of the system are. 

6. For the Interactive Graphics Version, the parameters of each 
feature mentioned above are available for modification as desired on 
an interactive basis. For example one may select specific values of 
various parameters and observe their effect on the resultant solution. 
The program is interactive and prompts the user in the choice of 
responses indicating his desires. 

There are two versions of this program. One version is called 
the Batch Version and runs on the IBM 360/67. The other is an 
Interactive Graphics Version and is implemented on the Electrical 
Engineering Department's XDS 9300 with the associated ADAGE 
Graphics Terminals. The interactive graphics version was primarily 
practical due to the power of the "intelligent'' Adage Terminals. 

In Chapter Ila general description of the program and its various 
parts are discussed. A major portion of the coding is the compiler 
and interpreter which are one of the primary features of this program. 

Chapter IV concerns itself with a description of the program from 
the users viewpoint. In addition to the general comments, there are 
included sections which deal separately with the Batch Version and 


the Interactive Graphics Version. 


el 





The use of the Interactive Graphics Version requires a knowledge 
of how to use the graphics system. This is explained by a suggested 
first session exercise. In this exercise the user is led step-by-step 
through all the techniques required to use the program. 

Chapter V ties all the above together by presenting a large number 
of actual problems that have been studied with this program. The 
input deck is explicitly presented along with the resultant output. 

A special benefit of the analysis feature of this program is its 
usefulness as a teaching or laboratory tool. This is especially true 
for the Interactive Graphics Version used to analyse a second order 
time independent system. For example the user could type in the 
description for a simple second order pendulum system. The slope 
lines give a very dramatic overall view of the effect of various types 
and amounts of damping. Since the system is interactive, the user 


can immediately pursue and verify theoretical observations. 


Z 





_u. GENERAL DESCRIPTION OF THE PROGRAM 


A. INTRODUCTION 

The program consists of five primary parts. The first four are 
the compiler, interpreter, slope line generator (used in the second 
order time independent case only), and the solution generator. Each 
of these will be explained in more detail later in-this chapter. For the 
present, a brief functional description will be given. 

The ''Compiler' takes the differential equation or system of dif- 
ferential equations, as the case may be, and through a process to be 
described later comes up witha polish string [Forsythe, 1969, Bur- 
roughs, 1964] of integer codes stored as a single vector. This polish 
string is the compiled version of the input system of equations. 

The 'Interpreter,'' when directed, uses the polish string and the 
current value of all its required arguments and produces a value for 
the derivative of each state of the system with respect to time. The 
required arguments are indicated by the system equations and may 
include the current values of the states, time, and the values assigned 
the constant coefficients A through H. 

The resultant state variable derivatives as defined by the inter- 
preter are needed by both the slope line generator and the solution 
generator. The slope line generator calculates the grid of slope lines 
as shown in Chapter V. The solution generator, as the name implies, 
generates the time solution. 

The fifth primary part of the program is the coordinator. It is 


the coordinator that ties the above together and makes it work. It 
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provides the required scheduling based on the results of the user's 
modifications. These modifications are accomplished on an interactive 
basis when using the XDS/AGT Interactive Graphics Version and via 


data cards on the IBM 360/67 Batch Version. 


B. MATHEMATICAL BASIS OF THE PROGRAM 


1. The General Mathematical Concept 


Assume an Fe order differential equation of the form 


dx 
= + f(x) = g(x, u) (2-1) 


dt 


where f(x) is a function of the n-1 derivatives, and g(x,u) is a function 
of the inputs and the n-1 derivatives. This equation can be rewritten 
as n first order differential equations. Let us assume the following 


notation and substitutions 


x = Z(1) 

x! Ze Za Ome) 

x! =a = "2 DOT (2) 

xt! = Z(4) = ZDOT(3) (222) 
oe So7in) = ZDOTMHe?). 


Rewriting equation 2-1 using the substitutions and notation specified 


in Equations 2-2 gives 


Te 


ZDOT(1) 


ZDOT (2) Z(3) 


(2-3) 
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ZDOT(n-1) Z(n) 


ZDOT(n) 


Piez, u) iGey 


=% 


Each of the Z's will be called the state variables of the system. 
Note that for the first n-l state variables the derivative of the fee state 
variable of the system is equal to the itl state variable. The deriva- 
tive of the we state variable is equal to g(x, u) - f(x). 

The solution generator needs only these derivatives of the state 
variables of the system. The compiler and sme ie: use the above 
observations as the basis for the method of implementation. 

Assume that equation (2-1) has been presented to the system. 


The first action of the compiler is to recognize the term which 


dXx 
dt® 


x'' followed by n apostrophes. This 


is denoted on the data card as an"! 


notation is completely analagous to the use of dots to denote time de- 
rivatives. The notation will be explained in detail in Chapter III. 
Having recognized this term it now knows the order of the differential 
equation. Next the expression g(x, u) - f(x) is compiled. The compiled 
version will be used by the interpreter to evaluate the desired 
derivatives. 

The interpreter first evaluates the compiled version of g(x, u) - f(x) 
which becomes the value of the highest order state derivative. The 
other derivatives are defined using the state values as defined in 
iq. iz=3 ). 

z2. Mathematical Derivation of the Slope Lines 

The slope lines are used only for second order time indepen- 
dent systems. Thus the comments of this section will apply toa 


second order system only. The general system can be written as 


x'' + f(x) = p(x) (2-4) 
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where f(x) + g(x,u) are as defined in equation 2-1. 
Letting x = Z(1) and x' = Z(2) gives 


tedro= AbDOT(1y= Z(2) 


dx'/dt = ZDOT(Z) = g(x, u) - f(x). = 
For a given time t the slope of the solution trajectory 
is 
slope = (dx/dt)/(dx'/dt) (2-6) 


ZWOT(2)/Z(2). 
Thus knowing a position and a corresponding slope the desired lines 


can be generated using simple trigonometry. 


C. COORDINATOR 

Figure 2-1 is a simplified block diagram of the program and out- 
lines the function of the coordinator. The coordinator provides the 
overall coordination and interface for obtaining the desired chain of 
events. The program is best summarized by reviewing this block 
diagram. 

The first action of the program is to obtain the system and plot 
description data. This includes the system equations, the plot para- 
meters suchas scales and sizes, and the solution information suchas 
the initial conditions, initial and final times. With this informationa 
complete run through the program can be made. 

Next the compiler compiles the system description data. The 
results of this compilation, if successful, are saved ina vector (a 
1-dimensional array or stack). If the result of the compilation or one 
of the steps which uses the compiled version encounters an illegal 


situation then an appropriate message is presented and the program 
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goes on to get new data. For the Batch Version the next new problem 
would be started. For the Interactive Version the user would make 
the ——_ corrections and proceed. 

After a successful compilation the plot is started and the grid of 
slope lines is computed and plotted. The time solution is then com- 
puted and plotted. 

At this time, the program is ready for ae next set of data from 
the user. Depending on what new data is obtained several courses of 
action are taken. If the system description equations are changed 
then a recompilation is required. The program flow after recompila- 
tion is the same as before. If the system equations were not changed, 
then the coordinator checks for a change in the constant coefficients 
or the plot parameters. A change in these will cause a transfer to 
the slope generating segment. 

The Batch Version differs slightly at this point. Instead of asking 
if the ''coefficients or plot parameters changed? ,'' the question should 
read: ‘''Is the next solution desired on a new graph?'' This allows a 
family of solutions to be plotted on the same graph if so desired. 


The above cycle of events continue until the user is done. 


De peo CRIPTION OF THE COMPILER 

The programming details of the compiler are not of central im- 
portance to this thesis. A copy of the coding is in the overall program 
printout included as Programs A and B. Frequent use of comment 
statements including explanation of codes, key variables, and such 
are included. 


Referring to the simplified block diagram, fig. 2-2, a very brief 


description of the compiler will be given here. The compiler accepts 
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the input a line ata time. The recognition of the elements of the 
language is accomplished by brute force. For example, let us assume 
that a "T' has been recognized. This T could be the variable T or 
the first letter of either ''TAN" or "THEN". A check of the following 
characters quickly reveals which of the three it is. Whena language 
element is recognized a corresponding integer code is placed ina 
stack. At the end of the recognition phase, the entire set of equations 
will have been coded into a string of integer codes stored in the stack. 

A blank card or line signals that the recognition phase has been 
completed. At this time the string is in the form known as "infix.''! 
The final task of the compiler is to convert this infix string toa 
"polish" string. The polish string has the advantage that it is quite 
easy to implement an interpreter for it. This interpreter will be 
discussed in the next section. 

A part of the compilation work is to recognize and use information 
such as the order of the differential equation. The order of the dif- 


ferential equation is determined from the equations and is required by 


the interpreter. 


Peet eeRIPIION OF THE INTERPRETER 

The purpose of the Interpreter is to take the compiled code which 
was generated in the compiler and values for the required operands 
which were specified by the user and to provide the derivatives that 
are needed by the slope and solution generators. The use of these 
derivatives will be explained in the next two sections. 

The functional block diagram is presented in figure 2-3. Upon 
entry to the Interpreter all the required operands have been defined. 


These are the present values of each of the state variables, the values 
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of the constants (A through H), and the problem time T. The Inter- 
preter takes one by one each coded entry in the polish stack (result 

of the compilation phase) and carries out the required operation. The 
various types of actions are summarized by the block diagram. 

The code may be carrying information such as the order of the 
differential equation and whether the differential equation to be used 
is a function of time. These codes cause the proper variables and 
Meresetowbe set. 

The logical branching is done in the following manner. Each 
statement in the stack is separated by a code representinga ";". If 
the interpreter evaluates the argument of an IF statement to be false, 
then it skips through the code to the entry following the next ";"'. 

Thus it has skipped the part of the statement indicating action to be 
accomplished if the logical argument is true. The SKIP instruction 
uses the same method but will skip several instructions as requested 
by the programmer. 

The arithmetic evaluation of the polish stack is accomplished by 
the use of software LIFO (last in first out) stacks. The basic prin- 
ciples involved in evaluating a polish stack are well known [Rosen 1967, 
Forsythe 1969]. For ease in implementation three different LIFO 
stacks have been used. One is used for the arithmetic operations, one 
for logical operations, and one to hold addresses. Thus they are 
respectively real, logical, and integer stacks. Associated with these 
stacks there are the required pointers to keep track of their top 
positions. 

When the interpreter notes that a differential equation has been 


evaluated it terminates the interpretation phase. As will be explained 
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in Chapters ]II and IV the program uses the first differential equation 
whose conditions for use are met. There it will be shown that the 
program easily handles a system whose behavior must be described 

by different differential equations in different portions of the solution 
space. The last action of the Interpreter is to define the desired 
system derivatives as described in Section B. These derivatives plus 
the order of the system are passed back to the calling slope or solution 


generator. 


a Peto CRIiPTION OF THE SLOPE GENERATOR 

The Slope Generator presents a grid of slope markers over the 
area of the phase plane. The use of these slopes is demonstrated in 
Chapter V. The mathematical basis of the slopes was described in 
Section B of this chapter. 

For each grid point on the phase plane the interpreter is called. 

If the system, for that grid point, is not second order or time indepen- 
dent then it skips onto the next grid point. The effect of this skip is 

to suppress the slope markers at that point. If the system is of second 
order and not time dependent then the slope is defined and is drawn 

at the current grid point. 

The resultant derivative information from the interpreter and the 
position information are used to calculate the desired slope. The 
slope is used in turn to compute the vector that represents the slope 
at that point. This derivation is made somewhat more complicated 
by the scaling options available to the user. Thus both the slopes and 
the vector length calculation must consider the horizontal and vertical 


scale factors. 
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Tew eoCRLPTION OF THE SOLUTION GENERATOR 

.The Solution Generator incorporates a very straightforward use 
of the Runge-Kutta method. The heart of the Runge-Kutta solution 
is the function subroutine RKLDEQ. This function subroutine is part 
of the IBM/360 library. For flexibility the actual cards are included 
in this program. 

After each solution point is generated the result is stored for 
future plots. Inaddition, the interactive graphics version presents 
the solution as it takes place. The solution rate for a typical system 
is about ten iterations per second with graphics. Thus the resultant 


presentation is approximately real time for a solution time step of 


0.1 seconds. 
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~ Ill, GENERAL DESCRIPTION OF THE LANGUAGE 


fe TNT RODUCTION TO THE LANGUAGE 

The language primarily resembles the familiar (to most) form of 
FORTRAN. Anyone familiar with one of the languages; FORTRAN, 
Algol, Basic, or PL/I will find the language used here quite easy. 
The requirements of this language are very limited and of a special 
nature. Thus, since we have a special purpose language, it is not 


complicated by the generality required of most languages. 


B. ELEMENTS OF THE LANGUAGE 
1. The Character Set 
The following characters are used in this language. 
a. The alphabetic characters as applicable. 
b. The numeric characters, 0 through 9. 


c. The special characters: 


x 


we 


mo -~_—™-? 


blanks 
d. Blanks are optional and may be used as desired for 
clarity and neatness except in the middle of a word, variable name, 
or number. The semicolon is used as a delimiter allowing more than 


one statement per line and may optionally be added at the end of a line. 
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2. Constants 

A constant is a fixed, unvarying quantity. The constants used 
in this program are of two general types. They can either be self- 
defining numeric values appearing in the source statements, or con- 
stants defined externally to the program. 

There are eight externally defined constants with names A, 
B, C, ... H. The values of these real constants are defined by the 
user externally to the system statements. Thus the values of these 
constants may be changed as desired, between executions, without 
recompilation of the systems equations. The method of defining these 
constants is explained in Chapter IV. 

The self-defined constants are recognized by the compiler 
during compilation. With the exception of integer powers all self- 
defined constants are treated internally as real. However, when 
used in the program, the use of a decimal point is optional. That is, 
the decimal point may be either explicit or implicit. 

In the case of exponentiation, an integer, positive or negative, 
used as a power will yield the expected integer power results. On the 
other hand, non-integer powers are computed using the natural log- 
arithm/exponentiation method. 

Examples: 

Valid real numbers: 9.05, 0.06, 78, etc. 
Valid use of integers: 9.7%**4, V12**-3, V1Z ** - 3, etc. 
Valid use of non-integer exponents: V23**0. 897, etc. 
3. Variables 
A variable is a named quantity whose value may change 


during execution of the program. Variables are specified by name. 
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The name must start with the letter ''V'"' followed by one or two 
decimal digits representing the number 1 to 40 inclusive. These 
variables are always considered as real variables. 

Example: VI, V23, V39, etc. 

In addition the letter ''T'' is used asa pe T is the real 
variable which represents the problem time; T is the time variable of 
the system solution. ''X'' is the variable of the differential equation 
and will be explained later. 

4, Expressions 

An expression is a string consisting of constants, variables, 
and function calls interspersed with operators and parentheses. The 
formulation of an expressionis governed by normal mathematical 
convention and the rules given below. There are two kinds of expres- 
sions; arithmetic and logical. The value of anarithmetic expression 
is a real number. The value ofa logical expression is always a truth 
value, TRUE or FALSE. Expressions may appear in If Statements, 
Assignment Statements and in Differential Equation Statements. 

a. Arithmetic Expressions 


(1) The arithmetic operators are as follows: 


Paikometic Operator Definition 
moe Exponentiation 
* Multiply 
/ Division 
+ Addition 
Subtraction 


(2) Rules for the Construction of Arithmetic Expressions: 
(a) All desired computations must be specified 
explicitly. For example, to multiply the variables Vland V3 one must 


write V1l*V2 and not V1V2. 
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(b) No two numeric operators may appear in sequence. 
(Note: ** ig considered as a single operator). Thus A* - V19 is 
illegal but may be corrected by the addition of parenthesis obtaining 
A*(-V19). 

(c) Parenthesis may be used as desired for clarification 
and to force the proper order of computations. 

(d) A mathematical function, together with its argument, 
may be used and is treated as a variable. 


(3) Order of computation in Arithmetic Expressions: 


Operation Hierarchy 
Evaluation of functions lst 
Exponentiation 2nd 
Multiplication and division ged. 
Addition and subtraction 4th 


(4) Examples of Arithmetic Expressions: 
(A*V3 + 5.7)/EXP(V1*B*T ) 
-A*(1,0 - X**2)*X!' - BeX 


A*X'*ABS(X') + SIGN(X) 


b. Logical Expressions 
(1) Relational Operators. There are six relational operators. 
They, and the logical operators which follow, must always start and 
end witha period. The relational operators are binary operators 
which operate on arithmetic expressions and express an arithmetic 
condition which can be either true or false. The relational operators 


are as follows: 


Relational Operator Definition 
oat. Greater than 
. GE, Greater than or equal to 
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ol: Less than 


we. Less than or equal to 
= . EQ. Equal to 
. NE. Not equal to 
Examples 
Vie. LT. xX! 


e742 .GE. EXP(A** - 0, 5) 
eee. “we. ORS X .GT. A 
(2) Logical Operators. The three logical operators are 
-NOT. , .AND. , and .OR. . Their use is defined in the 
ordinary sense of logic. 
(3) Rules for Construction of Logical Expressions: 

(a) Relational operators may only operate on expres- 
sions which, when evaluated, have a numeric value. 

(b) Logical operators may only operate on expressions 
which, when evaluated, have a logical value (i.e., TRUE or FALSE). 

(c) Relational operators may not appear in sequence. 

(d) Two logical operators may appear in sequence only 
if the second one is the logical operator .NOT. 

(e) A relational and a logical operator may not appear 
side by side in an expression. 

(f) As in Arithmetic Expressions, parenthesis may be 
added as desired for clarification and to force the proper order of 
calculations. 

(g) If the logical operator .NOT. operates on an ex- 
pression with more than one term, the expression must be enclosed 


with parenthesis. 
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(4) Order of Computations in Logical Expressions. 


Operation Hierarchy 
Evaluation of functions Ist 
Exponentiation 2nd 
Multiplication and division oma 
Addition and subtraction 4th 
ee eee e., .NE., .GT., .GE. Stak 
. NOT. 6th 
AND. 7th 
ne. 8th 


If two operators have the same hierarchy, then they are used in the 
order they appear left to right. 

(5) Examples of Logical Expressions. 
feces GCE, A) .AND. (X' .LT. 0) 
SecGen. 1.5 
mepeedee 157 OR, ({X .NE. 0.5) .AND. .NOT. (X .LY. X'") 
woe. 6 cle: COR, 6 6(XhlCUCUNE. «20.5)) . AND. . NOT. (X .LT. xX" 
These last two are not equivalent since the added set of parenthesis 
has altered the computation order in the last case. 

5. Mathematical Functions 

The language provides for the use of various mathematical 
functions. The present library is listed below. One calls these 
functions in an implicit manner by simply using the proper function 
name, followed by its argument in parenthesis, as an ordinary 
variable. The argument itself may be any valid mathematical ex- 
pression which in turn may contain function calls. 

a. Examples of Use. 


EXP(V1) 
EXP(V1*V3/SQRT(SIN(X/X'))) 
(SIN(X) + COS(X **2))/3. 75 
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b. The various available functions and their preassigned 


names are expressed below. 


Name Description of Function 

SIN Sine of an angle in radians 

COS Cosine of an angle in radians 

TAN Tangent of an angle in radians 

ABS The absolute value of the argument 
EXP Exponential of the argument 

LN The natural logarithm of the argument | 


LOG The base 10 logarithm of the argument 
INT Truncates the real number to an integer then reconverts 
Mm toa real. 


SIGN Magnitude of 1. 0 with the sign of the argument. 
6. Arithmetic Assignment Statements 
The goal in developing the above language elements was to 
make use of them as the building blocks for Arithmetic Assignment 
Statements, Control Statements, and Differential Equation Statements. 
This section will discuss the Arithmetic Assignment Statement. 

The end result of using an Arithmetic Assignment Statement 
is to compute a new value fora variable. The general form of an 
Arithmetic Assignment Statement is a = b 
where: ais a valid variable name (written without a sign). 

b is any mathematical expression following the above 

guidelines. 

= is an order to compute the value of the expression b 

on the right, and then give that value to the variable 
named on the left. 
A lengthy discussion can be made onthe difference between the 


Arithmetic Assignment Statement and a mathematical statement of 
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equality. It is perhaps unfortunate that an '"'="' is traditionally used in 
both.cases. However if one keeps the above definition in mind, there 
should be no confusion. For someone completely new to computer 
programming methods, it is recommended that one of the many good 
fundamental FORTRAN manuals be consulted. The formation of and 


use of the Arithmetic Assignment Statement is identical with that of 


FORTRAN. 
Examples: 
] a 
let Vl =— + —— 
2x 2 
4x 


then the coded version is: 


V1 = -1. /(2. *X) + A**2 /(4*X 2) 


3 2/5 
Let V5 = (A + Bx ) 


then the coded version is: 
V5 = (A + B¥X%**3)%#*(2, 0/5. 0) 
i 
Let V9 = (SIN“x + cos’x)? 


then the coded version would be: 


Me =e(SIN(X)**2etCOS(X)**2 )**0, 5. 


tan * | 





Inetey 25 fee + log 
cosx 


then the coded version would be: 


V23 = 1. 0/COS(X) + LOG(ABS(TAN(X/2. 0))). 


Let V5 aS. + Be 
new old 


then a coded version would be: 


V5 = V5 + ABS(A*X*X). 
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7. Control Statements 
a. The Skip Statement 
The control statements are of two general types. One 
is a simple unconditional transfer statement called a Skip Statement. 
Its form is: 
SiclP +N 
where: N is an integer that specifies how many statements should be 
skipped. 
b. The Arithmetic IF Statement 
The other more powerful control statement is the 
Arithmetic IF Statement. This is essentially the statement that one 
is familiar within FORTRAN. Its form is: 
IF(a) THEN b; 
where: ais any valid logical expression as previously defined. 
b may be an Arithmetic Assignment Statement, a Skip 
Statement, or a Differential Equation Statement. 

It commands the computer to evaluate the logical expressiona. Ifa 
is true then do b, otherwise skip b and go on to the next statement. 
The use of the keyword THEN is optional as is the "';". 

c. Examples of Valid Control Statements 


SKIP 3 

ime .Gr. xX" ) THEN SKIP 4 

mix LE. ABS(X')) THEN Vi = A*xX*%2 

Pee GE 0.0)) THEN X'' + 0.2*X' + X = 0.0 


8. Differential Equation Statements 


a. Notation 
Differential Equation Statements are intended to be as 


natural as possible. The apostrophe notation, which is widely used, 
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has been adopted. Thus this program used x'!' as the symbol for the 
second derivative of x. 

| Other derivatives of X are handled in an analogous manner. 
With the exception of this notation for the derivatives, the remainder 
of the differential equation is written using the SonsteNenen methods 
used for mathematical expressions. Note that in the case of a Dif- 
ferential Equation Statement, the statement actually expresses an 
equation. Thus for this statement only, the equal sign denotes 
equality rather than an assignment. 
Pemaple: X'' + A*X'*A BS(X') + X**2 = VIST 

b. Writing a Differential Equation Statement 
The form of a Differential Equation Statement is: 


aye 


a 
where: a denotes the highest derivative of X written using 
the notation of the previous section. 
b and c denote any valid mathematical expressions 
using any or all of the constants and variables and 
using any or all of the lower derivatives of X as 
variables. 
It can be seen that the rules for writing b and c are extremely liberal. 
Thus the type of differential equation that can be handled is very 
general. The differential equations may be extremely non-linear, 
and have coefficients that vary with time and/or position. 
ec. Examples 


The above generality can best be demonstrated by a few 


examples. Valid constructions are: 


Peet ee )~CU tT COL 75 *X' 6+ «(20.0 
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Toe + (OG )=eX + SIN(X) =.10..0 


ei ttt + Kt 4 XU/xX = LOG(X!!'! #0, 5) 


x 


Cc. USING THE LANGUAGE 
1. General Comments 

The goal of all the above elements of the language has been to 
develop tools for use in defining the differential equations of a system. 
The above language features make it possible to handle a wide variety 
of system types. This section demonstrates, by use of several ex- 
amples, how to use the elements of the language to describe various 
systems. The series of examples will start with simple applications 
and then progress into more involved descriptions. 

2. A Simple Differential Equation 

The simplest, but still very useful and powerful, form of 
system description consists of one Differential Equation Statement. 
Let us take for example Van der Pol's 
Equation: X'' - A*(1,0 - X**2)*X' + B¥X = 0.0 
This one data card, or line on the graphics terminal, completely 
describes the differential equation with the exception of A and B. 
The variables A and B are specified on an additional data card. When 
it is desired to alter the values of A and B they are changed independ- 
ently of the system equations. Considerable execution time is saved 
by eliminating recompilation of the system equations when it is only 
desired to alter the value of a coefficient. If, however, one did not 
anticipate varying coefficients A and B, then the differential equation 
could be coded as follows, for example: 


pee 0, 95(1, 0 - Kek2)eX' + 1.7*X = 0.0 
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Assume now that it was desired that the first coefficient be 0.90 vice 
0. 93. There is still no problem, the required recompilation is trans- 
parent to the user. The on line user will only notice a somewhat 
longer delay time before the solution starts. 
3. More than One Differential Equation Statement 

In most control system applications more than one differential 
equation is needed to describe a given system. However, only one 
differential equation is applicable at a given moment in the solution. 
For example, due to a particular type of non-linearity, the selection 
of the proper differential equation may be dependent on the sign of the 
velocity or position. In another situation, the selection of the dif- 
ferential equation may be a function of time. Numerous examples of 
various control system descriptions appear in Chapter V. 

For purposes of illustration, we shall assume the following 
system: 


If X is less than -0.2 
then the differential equation is 
ee OZ et OX =. 0.2 
If X is greater than 0.2 
then the differential equation is 
re tt OS 2x UF UX C= -0.2 
If X is neither of the above 


then the differential equation is 


met 0.2 Ut CX C= CO 
The above system could be coded as 


mec, - 0.2) THEN X'' + 0.2*xX' + X =0.2 
mio eeT. 0.2) THEN X'' + 0.2*xX' + X = -0.2 
Paes tt Oeerx’ + X = 0.0 
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During the execution of the solution, when the differential 
equation is referred to, the Interpreter will test the argument of each 
IF statement starting at the top. It will use the differential equation 
corresponding to the first condition it finds to be true. Note that the 
third statement is not preceded by a condition. The effect of the third 
statement is that the condition for its use is always true. Thus if the 
Interpreter does not use one of the first two differential equations it 
will use the third. 

To further clarify the selection of the right differential 
equation, consider the following: 
ree orzstx' + X = 0.0 
me TD, =0.2) THEN X'' +_0.2*X' + X = 0.2 
ie .Gl. 0.2) THEN X" + 0.2*xX' + X = -0.2 

This is the same set of statements as before except that the 
bottom statement has been moved to the top. However, the result is 
drastically different. The condition for use of the present first dif- 
ferential equation is true by default. Thus it will be used regardless 
of the value of X. Again, the point to be stressed is, only one dif- 
ferential equation can be used at a given moment; the one used will be 
the first one witha true condition. 

4. Using the Remaining Features of the Language 
Often it is convenient to be able to define and calculate some 
variables before the Differential Equation Statement. Consider the 


following coded system description: 


IF((4. 0#X + X') .GE. 0.0) THEN X'' + 0.02*X' + 0.3 
irene x) | LT. 0.0) THEN X'' + 0. 02x" - 0.3 


1] 


Using features of the language previously described this could be 


recoded as follows: 
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femeet 0.5, IP(arorx + X") .LT, 0.0) THEN VI = -Vl 
me + 0.027X' + Vi = 0.0 

The semicolon allows more than one statement per line. 
This recoded version has the same effect as the previous version. 
Note that there is only one differential equation. The V1 used in the 
differential equation is defined in the line above in such a way that 
the effect is the same as the two differential equations of the first 
version. 

For this system an even simpler coding is possible. 
X'' + 0.02*X' + 0. 3*SIGN(4. OFX + X') = 0.0. This example 
demonstrates the versatility of the language by showing the same 


system coded in three different ways. 
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IV? -PRNCIPES AND METHODS @P USE 


A, GENERAL COMMENTS 

The previous chapters have concentrated on various specific 
aspects of the program. In Chapter IJ the internal structure of the 
program was described. Chapter III described the language associated 
with use of the program. This chapter will endeavor to pull these to- 
gether with specific instructions on the use of the program. In Section 
B those areas common to both the Batch Version and the Interactive 
Graphics Version will be discussed. Sections C and D will deal with 


the peculiarities of the two versions respectively. 


B. GENERAL PRINCIPLES 
1. System Description 
One of the first tasks for the user is to encode the system 
equations in the form usable by the computer. One of the primary 
features of the program is that this is done ina very natural way. 
For simple systems this will involve writing only a single equation. 
For systems with a complicated control strategy the user may have 
to write what amounts to a small program to completely specify it. 
This general area was covered with examples in Section C of Chapter 
lil. Chapter V will further amplify this with a large 
number of complete examples. 
2. Use of the Constants 
The use of the constants was described in Chapter Hl. How- 
ever for the completeness of this Chapter they will be mentioned 


again. Often when preparing to study a particular system it may be 
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desired to vary a parameter over several values. One way to accom- 
plish this is to vary the numerical value in the system equations. 
However this amounts to starting a whole new problem and the system 
equations must be recompiled. 

An easier way to accomplish this parameter variation is to 
use the constants. These are the constants labeled A - H. When 
writing the system equation use the constants in place of an actual 
numerical value and define them separately as desired. When it is 
desired to change the parameter, redefine the associated constants. 

Another handy use of these constants is when the same 
number is to be used several places in the system description. In 
this use a constant would be defined once and then used as often as 
desired. 

3. The Phase Plane Window 

The output of the program is a phase plane. The slopes and 
solutions are presented on this phase plane. This phase plane can be 
visualized as a window in which the user has control over the place- 
ment, size, and scales to be used in its construction. The controls 
are independent in the horizontal and vertical directions. Following 
are the control parameters and their use: 

X-SCALE The horizontal scale in units per inch. 

Y-SCALE The vertical scale in units per inch. 

X-CENTER The X coordinate of the center of the window. 

Y-CENTER The Y coordinate of the center of the window. 

X-SIZE The horizontal dimension in inches. 


Y-SIZE The vertical dimension in inches. 
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4. Solution Parameters 

The remaining values to be specified are the parameters of 
the solution. These include the initial time, time step, final time, 
and the initial conditions. 

In most cases, the user will let the initial time be zero and 
the final time be some positive value with an associated positive time 
step. However, this is not required. In fact the final time may be 
less than the initial time. Inthis later case the time step would be 
negative. 

The program can handle systems up to eighth order. Thus 
there are provisions for specifying up to eight initial conditions. 
This limitation to eighth order was arbitrarily set by the ease of 
handling the initial values on input. Eight initial values fit nicely on 


one input card. 


eee oesCirPic DETAILS ON THE BATCH VERSION 

1, Planning for the Batch Run 

For the batch version the user must plan the entire run 

beforehand. One may be uncertain of some parameters of the 
problem onthe first run. If this is the case, the user should take 
advantage of the flexibility of this program on the first submission. 
That is, the first submission should include combinations of values 
such that when the results are returned one will know what values 
to choose for the next run. For example, take the problem of the 
choice of time step and final time. Here it would be good to make 
a guess at a suitable time step and final time. Then run an addition- 


al two solutions with these times multiplied by ten and divided by 
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ten. The added effort for these additional solutions would be trivial. 


When the results come 


the best values. 


back it will then be quite easy to decide on 


The same idea canapply to other areas such as scaling. The 


time required to punch the added cards for a few extra runs can save 


considerable time lost to turn-around-time if an added exploratory run 


is required. 


oe Data Card Arrangement 


Cards 1-2 


BOR; Os .. 


Plot title including your name (6A8) 
Format 


X-SCALE(Units /inch), 

X-CENTER (Coordinate of X-CENTER) 
X-SIZE (Inches) 

Y-SCALE (As above) 

Y-CENTER (As above) 

Y-SIZE (As above) 

These six items are placed on one 
card in the above order using a 

(6F10. 0) Format. 

if (X=SIZE .GT. 9.0) thenmmcne ploris 
rotated CCW onthe paper by a quarter 
turn. Due to paper size only one of 
the dimensions may exceed 9. 0 inches 
and must not be more than 15. 0 inches. 


SIZE of slope markers (inches, 0.2 
works nicely), NUMBER per inch (3 
recommended) 

In (F10.4, 11) Format 


BOUATIONS TO BE USED: Following 
the rules of Chapter III (80Al Format). 


Then follow four cards per solution desired as 
described below: 


-First- 


-Second- 


-Third- 


VALUES of A, B, C, D, E, F, G, H (8F 10. 4) 
Format 

Integration info (INITIAL TIME, TIME 
STEP, and FINAL TIME) (900 Solution 
steps Max), (3F10.4) Format. 

INITIAL CONDITIONS for Xe eee, 

.. . in(8F10.4) Format 
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-Fourth- A "WHATS NEXT CARD" use the following 


code: 
QO - There are no requests following. 
1 - The following are the second and 


third cards with new integration 
and initial condition information 
only. I want this new solution to 
be plotted on the present graph. 

2 - The following is another four cards 
with data for another solution of the 
present system on a new graph. 

3 - The following card starts a 
completely new problem starting 
with the plot title cards. 


3. Sample Data Deck 

Following is a typical sample deck of input cards. The com- 
ments enclosed within parenthesis are not part of the cards but are 
included in this listing for clarification only. These cards accomplish 
the following. 

The system X!'' + A*X' + SIN(X) = 0 is compiled and two 
solutions are obtained with A = 0.3. Then A is changed to 0.7 and 
two more solutions are run using the same time information and 
initial conditions as with A = 0.3. Thena completely new problem 
(a saturated servo system) is defined. 


The input cards used are: 


NELSON, H. G. Pendulum (GRAPH TITLE - lst LINE) 
oo = Aa! = SIN(X) = 0 (GRAPH TITLE - 2nd LINE) 
2.0 0. 0 8.0 2.0 0. 0 6.0 (WINDOW INFO) 
Oe2 4 (SLOPE MARKER SIZE AND NO. ) 
X'' + A*®X' + SIN(X) = 0 (SYSTEM EQUATION) 
(END OF SYSTEM EQUATIONS) 
0.3 (VALUE OF A) 
0.0 0. 05 2570 (TIME INFO) 
=e 0 a0 (INITIAL COMP Ions, 
] (WHAT'S NEXT SIGNAL) 
0.0 0.05 25.0 (TIME INFO) 
=O 5. 0 (INITIAL CONDITIONS) 
2 (WHAT'S NEXT SIGNAL) 
0.7 (VALUE OF A) 
0.0 0. 05 25.0 (TIME INFO) 
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7, 0 3.0 (INITIAL CONDITIONS) 


1 (WHAT'S NEXT SIGNAL) 

0.0 0.05 25.0 (TIME INFO) 

Ha) 5.0 (INITIAL CONDITIONS) 

3 : (WHAT's NEXT SIGNAL) 

NELSON, H. G. 0902 (GRAPH TITLE - Ist LINE) 
SATURATED SERVO SYSTEM(GRAPH TITLE ~ 2m LINE) 

0.3 0.0 6.0 0,2 0.0 8.0 (WINDOW INFO) 

O.2 4 (SLOPE MARKER SIZE AND NO, ) 


nec t T. (1) 0=2. 0)20.2 THEN X'' + 0.2*X' + X = 0,2 
IF(X, GT. 0.2) THEN X"! + 0.2%X' + X = -0.2 (SYSTEM EQUATIONS) 
Pat 0. 24K! + X = 0.0 
(END OF SYSTEM EQUATIONS) 
0 0.05 30. 0 (TIME INFO) 
7 (INITIAL CONDITIONS) 
(WHAT'S NEXT = BLANK = 0, THUS 
END OF DATA) 


QO. 
0. 


wee oe GCIFIG PE EAILS ON THE INTERACTIVE GRAPHICS VERSION 
1. Planning for the Interactive Session 
Since the Interactive Version is very easy to use and, as the 
name implies, is interactive the required preparation is minimal. 
This statement assumes the first-time user has read this manual. If 
a mistake is made or a bad guess is made for a parameter value, it 
can normally be corrected before the run starts. If one notices the 
error after the run has started, there is stillno problem. The pro- 
gram has a large amount of protection built into it. Hopefully, no 
mistakes or errors by the user can cause the interaction to terminate. 
Regardless of the immediate results of a mistake, control will return 
to the user for desired corrections and/or modifications. 
2. Check-off List for Loading the Program 
Following is a step by step list to be used to load and run the 
program. 
a. Turn onthe XDS 9300 Computer (if not already on) using 


the instructions attached to the console. Also turn on Sense Switch 2Z. 
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b. Mount the binary tape of the program on one of the tape 
drives and advance it to the load point. Set the "unit select'' to 0, the 
"density select" to 556 and the mode switch to automatic. The "unit 
power", "unit ready", "file protect", and "load point" lights should be 
on. 

c. Place the small card deck as specified in Appendix B in 
the card reader and press the "power on" and ''start'' buttons in that 
Gader. 

d. Ready the printer. Onthe XDS 9300, verify that ''Sense 
Switch 2"' is on. Then press the following buttons: 

HALT 

Rees ET 

CLEAR and CLEAR FLAGS simultaneously 
RUN 

CARDS. 

e. The card reader should begin reading cards and then the 
Tape Drive will begin operating. It will take about three minutes to 
load and link the program. 

f. Following the instructions at the AGT's, turn on the AGT 
unit to be used and load ''Gated". 

g. Verify that ''Gated'' is ready by pressing the "Text Edit" 
button on the moveable button panel. The following should appear at 
the bottom of the screen: 

Jie eb iOGCK SELECT MODE Bleek. 1" 

h. When the XDS console typewriter types "ENTER 

IDEV = 1 or 2" do the following: 

1. Ensure the ACT is ready 
2. Note the number of the AGT 
3 


Type "IDEV = 1*"', return 
or "IDEV = 2%", return as applicable. 
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i. This completes the preparation phase. Although this list 
appears long or complicated it is quite easy once it has been done a 
few times. Quite often the user will find both the XDS and the AGT 
already turned on and ready to go. It is strongly recommended that 
the initial use of this program be made during working hours when 
advice and/or help is available. 

ee Witerecting with the Program 

The interaction starts with a question to the user. ''Do you 
desire to start a completely new problem?, answer 'YES' or 'NO'", 
Answer "YES" or "NO" on the keyboard and then hit the return key. 
If you answer ''NO"', a problem already in the computer will be dis- 
played. This is useful for becoming familiar with the program. It 
is suggested for the first time user that the answer should be ''NO" 
to that question. Then the new user should try various manipulations 
on the given problem until a working knowledge of the program is 
obtained. 

After each run the control returns again to this same question. 
After the first pass the answer to the question gives the user the 
option of making modifications to the present system already on the 
screen or starting a complete new problem. If there are only a few 
changes to be made it is best to selectively edit the present problem. 
If it is desired to start a new problem or make many changes, it will 
probably be faster to use the 'INPUT'' mode. When using the "INPUT" 
mode, the system prompts the user for the desired information. 

The easiest way to learn the system is to use it. Thus the 
following section is a step-by-step procedure for using the program 


on several specific problems. It is suggested that the user closely 
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follow this guide using the given problems. Once the user is familiar 
with the procedures for using the program and the graphics terminal 
then he can proceed on his own. 

4. Suggested First Session Example 

a. Load the program as described in Section 2. The question 
me YOU DESIRE TOSTART A COMPLETELY NEW PROBLEM?, 
ANSWER YES OR NO" will appear on the screen. 

b. Type ''NO", return. At this time the internally stored 
program will appear onthe screen. In addition will be the message 
"MAKE DESIRED CHANGES AND CHOOSE NEW INITIAL CONDITIONS 
(TEXT OR LIGHT PEN) TO INITIATE NEW RUN". The signal to the 
program that the user is ready to proceed is to edit the first line of 
the initial conditions or to select initial conditions graphically with 
light pen. 

c. If satisfied with the initial conditions already on the 
screen, make the following null edit. Press the following buttons in 
succession: 

TEXT EDIT 

NEXT BLOCK - Check to see that the second line characters are 
slanted. If not use the NEXT BLOCK and/or 
PREVIOUS BLOCK until line two is slanted. The 
slanting indicates that line two is ready to be 
edited. 

GO # DIT - Corrections can now be made as desired to line 
two. For this time a null correction will be 
made, i.e., no correction. 

TERMINATE EDIT 

The program now has all the information it needs to make a run and 

the above edit of line two has indicated that the program is ready to 
be run. When the run has terminated the question of step a. will be 


repeated. A complete loop in the program has now been completed. 


d. Again answer by typing ''NO", return. 
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e. Before this next run the goal will be to change the initial 
conditions to -3.8 and 3.0 respectively. To do this "'select" and ''go 
edit'' on line two as in step c above. This time a correction will be 
made. While holding the typewriter key marked ''CTRL", space over 
with the ''F"' (forward) key until the underscore is immediately to the 
right of the first initial condition to be corrected. Now release the 
"CTRL" key. Press the 'RUB OUT" key the desired number of times 
to eliminate the portion to be changed. Then type in the new first 
initial condition. Be sure to replace as many characters as were 
removed. This keeps the whole line in the proper alignment. Next, 
Space over and correct the second initial condition using the same keys 
as used to correct the first initial condition. When you are satisfied 
with the line press the terminate edit button or the "RETURN" key on 
the typewriter. Either one will terminate the edit and initiate the next 
run. The end of this run will bring you back to the question of step a. 

f. Again answer by typing ''NO", return. 

g. The goal of this run is to select a set of initial conditions 
Hemnesthe light peny Press in*sequence the ''graphics edit", "go edit’, 
and ''cursor on'' buttons. Insure that the light pen is turned on by 
checking for the pilot light on the box at the end of the light pen cable. 
Take the light pen, point it at the cursor, press the little white button 
on the light pen and lead the cursor around the screen. When you are 
familiar with positioning the cursor, select with the cursor a desired 
initial condition, then press the ''move"' and then ''draw" buttons. 

With the light pen move the cursor to another initial condition on the 
screen. There should be a dot remaining at the location of the first 


initial condition pair selected. Pressing the ''move", "draw" buttons 


48 





will select the present cursor position as a second initial condition 
pair. This can be repeated up to fifteen times if desired. When you 
are through selecting initial conditions, terminate the edit. This will, 
as when selecting the initial conditions numerically, initiate the run. 
A solution will now be obtained for each of the selected initial con- 
ditions and then control will return again to the question of step a. 

h. Again answer ''NO", 

i. This time change A to 3.0. To dothis press "text edit" 
and then press next block until the characters in the line containing 
the A tilt. You are now ready to edit this line using the method of 
step e above. 

j. When through editing A initiate the next run using the 
technique of steps c, e, or g above. 

k. Once more answer ''NO". 

l. This time it is desired to edit the system equation to 
obtain X'' + X*X' + X =0. Press the ''text edit'' and then the ''next 
block" button about eight times until the equation characters slant 
indicating they have been selected. Make the modification as before. 
You do not need to replace as many characters as you removed when 
editing the equations. When ready initiate the next run using one of 
the previously discussed methods. 

m. When the solution has been completed and the question 
of step a returns answer ''YES". In response to each question make 
responses as suggested below. After each response hit the return 


ery. 


a) 





SYSTEM EQUATIONS: 


Vl = 4. O*X + X' 

eG. Oe tigN XY + 0, O2*X' + 0.3 = 0.0 
ma’ + 0. 0O2*X' - 0.3 = 0.0 

blank line 


X-SCALE 


r 

—N 

C) 

> 

iS 

GJ 
non WoW oa 
CONC ON 


ee 
INITIAL TIME 
TIME STEP 
FINAL TIME 0. 
PIRST I. C. .6 
remainder of I. C.'s, hit return, 


it return, will default to zero 


oOo © 


= 
0 


Ooroonwomnoondco 


n. Now the program will give the question ''MAKE DESIRED 
CHANGES (TEXT OR LIGHT PEN) TO INITIATE NEW RUN". At 
this point you can proceed as in steps c. through l. above as desired. 
Thus if any corrections need to be made as a result of incorrect 
responses to the above questions they can be made at this time. When 
ready to initiate the run use the method of step c. It is best to select 
the initial conditions numerically for the first run after changing the 
window parameters. Initial conditions selected with the light pen take 
on the corresponding values of that point relative to the current 
picture on the screen. 

o. You have at this time exercised all the control functions 


needed for full interactive control of the program. 
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Veo eee Ee een ttONS [TO VARIOUS PRQOBLEMS 


A. GENERAL COMMENTS 

The purpose of this chapter is to demonstrate the use of the 
program by using it to study several classes of centrol system prob- 
lems. It is not the intention of this chapter to derive the applicable 
equations, but rather to show the use of the program in their solution. 
If the purpose of the system is to reduce the error to zeroora 
minimum, then it is convenient to use error as the dependent variable 
instead of the regular output variable. 

In the following examples the problem will be defined first then 
the data cards will be shown, followed by the resultant phase plane. 
Although these examples will have been produced using the hard plots 
of the Batch Version, the input data is identical with that of the 


Interactive Graphics Version. 
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B. VARIOUS SINGLE EQUATION EXAMPLES 
Example 1 
Van der Pol's Equation 


X'' - Ax*(1,0 - X**2)*X' + BEX = 0.0 


WHERE: 
A=1.0 
B= 1.0 


FLEOT PARAMETERS: 


ASCALE = 0 voCATE 
ACENTER = 0.0 YCENTER 
POZE = 6.0 YSZ 


Size of slopes = 0.2 


Number of slopes per inch = 4 


Pie INPUT CARDS USED: 


NELSON, H. G. 
pan tt SO | X42 )eK’ + BEX = 0.0 


bs 0 0.0 6.0 Peo 0. 0 8. 
OeZ 4 

Xt - Ax(1,0 - X**2)%X' + BX = 0.0 
br, 0 iO 
0. 0 0.1 20. 0 
0.01 Ce 
] 
0.0 0.1 20.0 
a7 0 0.0 
l 
Oo. 0 Om A. O 
Ze 5 30 
1 
0:0 Out 20. 0 
ZO -3.3 
] 
0.0 On! 20. 0 
-2. 5 a. 0 
] 
G20 0.1 20. 0 
2010 =2, 0 
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X'' - A*(1.0 - X*¥*2)*X' + B¥X = 0.0 


Example 1: 





Mose. = 1,0 units/inch 


Pee. 5-1 


Y-SCALE = 1.0 units/inch 
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+ ArTX*X" + BeX = 0.0 


WHERE: 
A=1.0 
B=1.0 


Example 2: 


Eeol PARAMETERS: 


XSCALE 


XCENTER 


MoIZE 


Size of slopes 


Number of slopes per inch 


ma INPUT CARDS USED: 


NELSON, H. G. 


Ore 


X'" + AXX2X' + B¥X = 


K'' + AxXX%*X!' + BxX = 0.0 


Pe 0 OO 

Bec 4 

1.0 0 

oO 0. 02 
Lv 

] 

ten) O02 
220 

] 

0.0 0. 02 
3.0 

] 

0.0 QO. 02 

Za.) -2. 0 

] 

0.0 0. 02 

2n5 =o. 9 


6210 


L520 


moe 0 


ieee O 


15.0 


15,0 


YSCALE 


MOCENTER 


Voz 
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meemaple 2; 20'' + AtKeX' + BYX = O. 0 


ew) 


B= 


iO, 


A 


pence ae eae -_ ae aioe Sal ae _ 


ee N 
i ro 
ee = ee ee ee. 2 eee et, 7 
Bien as eee i eae — i 








= — ee —_——_ 1 ee ee iene oe 





f 

ai ‘ 

/ 4 
4 
, 





—— 


een an 
a a 


a ae 


8 nie tee rare = 
lee ‘ 
ap ee 


c 








Mege NILE = 1. 0 units/inch 


5-2 


Fig. 


vec = 1,0 units/inch 
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Example 3: 


x + (1. 


PLOT PARAMETERS: 


XSCALE 


Oiae 


XCENTER 


XSIZE 


mre )rXx' + X = 0,0 


Size of slopes 


Number of slopes per inch = 4 


6. 


ae 


The input cards used: 


MmeLOON, H. G. 
X'' + (1.0 - X**2)"X' + X 
6.0 


r. 0 
Dee 


mi (1,0 - X¥22)*xX' + X = 0.0 


e 
Mm © ore" 


OM NOKFNOKK OF KOO 
© 


t 
N 
© 
in 


0.0 
4 


Om0m 


0.02 


0.01 


0.61 


0 


0 


0 


YSCALE 
YCENTER 


YoIZiE 
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Peeommpreucae ' +t (1.0 - X¥%2)*X' + X = 0.0 





ae 


Q 
oO 


SNS 
Sh Se 
Sa 


oe 


Wee fo 
ee ere 
ae 


mie 5-3 X-SCALE = 1. 0 units/inch 


Y-SCALE = 1.0 units/inch 


Sf) 





Example 4: 
het Xoree te xX = 0.0 


PeOT PARAMETERS: 


XSCALE = Zea YSCATE 
XCENTER = 0. 0 YCENTER 
XSIZE = &. 0 eS 


Size of slopes = 0.2 


Number of slopes per inch = 4 


mae INPUT CARDS USED: 


NELSON, H. G. 
X' + X'ee2 + KX = 


i 


Z. 0 
0.2 


0.0 
+ 


0.0 
8. 0 


beet X42 + X = 0.0 


ms = he 
; OO ' 76 BSFe eo 
oO oO oO 


Ww * 
=, 
) 


© 


OO 1 Oe OF igo KE oOorf! 2S o oe 
8 * ja ° 8 ® ® * 8 
. © 
© 


a) 


05 
Za. 5 


0.05 


nN © 


40.0 


40.0 


40.0 


40.0 


40.0 


40.0 


40.0 


l. 
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Et ee XC EZ X= 0 0 


Example 4 


§ 
pore ee ee, ee 


oe 
-+— — — —— ae 


aq. Corea 


i 








any 
Se 
_— 
ee, 
~——_, 
| 
ane. me. ea, 
t 
mi eee ee eT ——, 
| 
aN 
¢23 
ey — — 
a a mn are ae neem ace 
can ate ee ene tee me —~—— nan, 





a 


Fig. 5-4 X-SCALE = 2.0 units/inch 


Y-SCALE = 1.0 units/inch 





Eeceimple 5: 


X!! + ASX + BEX=#33 = 0.0 


WHERE: 
A = -1.0 
fee) 0.25 


PLOT PARAMETERS: 


XSCALE = 0 YSGALE = 1.0 
XACENTER = 0.0 YCENTER = 0.0 
XSIZE = 8. 0 Yoiz ke = 6.0 


Size of slopes = 0.2 
Number of slopes per inch = 4 
Ge INPUT CARDS USED: 


NELSON, H. G. 


X'l 4 AXX + BEX#83 = 0,0 
1.0 0.0 8.0 1.0 0.0 6.0 
0,2 4 

X' 4 AXX + BEX#3 = 0.0 
_1, (ES 
0.0 0. 02 18. 0 
1.0 ie 
1 
0.0 0. 02 18.0 
1.0 2.0 
1 
0.0 0. 04 30.0 
0.0 0. 03 
1 
0.0 0.02 18,0 
0,3 
1 
0.0 0. 02 18.0 
0.3 
1 
0.0 0. 02 18.0 
1.0 
1 
0.0 0. 02 18.0 
A 0 
1 
0.0 0. 02 18.0 
“1.5 
4.0 0. 02 18.0 
ins 
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eee ee 


ee ——— 
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— 
= _ eee eee 


yout/sjztun O°] 


yout/szIun O°] 





Be ~ a 
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: : ~ os -- a 
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Example 6: 
pee + (170 - ABS(X))*xX' + X = 0.0 


PROT PARAMETERS: 


XSCALE = boo oer: 
XCENTER = Ome YCENTER 
ASIZE = 6.0 Worse 


Size of slopes = 0.2 


Number of slopes per inch = 4 


tas INPUT CARDS USED: 


MIEESON, H. G. 
X'' + (1.0 - ABS(X))*X' + X = 0.0 
l. 


FeO OU 6. 0 0) O20 
0.2 4 
x £ (1.0 - ABS(X))*xX' + X = 0.0 
O70 
0. 0 O02 15. © 
1.0 
] 
0. 0 0. 02 LS. @ 
nD 
1 
0. 0 O02 Pe 
ZO 
) 
0. 0 OrOZ PS. 0 
2.0 
1 
Oc eee 15. @ 
-2.5 
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fe Oe ABS OX))*X' + X = 0.0 


Example 6: X"' 


eG ee 
SN VS oe 
ae 
|---| +++ 
rol os 
f | ht 
PP id dad 
ee 7 7 
1 jf jit 
Guy fae te 
Vy a 
Yi 
ee 

‘ee ae 

-— =f. 
SNK NS 
RN 
SS 


» - 
WS ~~, 


Ss 


oo 








ee 








Fig. 5-6 X-SCALE = 1.0 units/inch 


Y-S'ALE = 1.0 units/inch 
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Example 7: 
Ps oN TP Pero xX) = 0.0 


Peet PARAMETERS: 


XSCALE = Zao 
XCENTER = O70 
XSIZE = Orne 


Size of slopes = 0.2 


Number of slopes per inch = 4 


Dre INPUT CARDS USED: 

MeloON, H. G. 

i + X' + X*ABS(X) = 080 
Z, © 0.0 6.0 
0.2 = 

a ' +X’ + XTABS(X) = 0.0 


0.0 

0.0 0.01 8. 0 
a. & 

] 

0.0 Cup 500 
-3.0 

] 

0 Oot 8. 0 
0 

1 

0.0 0.01 8. 0 


Wn e 
mn 


mo CA LE 
YCENTER 


YolizZEe 
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= 0.0 


Example 7: X'! + X!' + X*¥ABS(X) 





2.0 units/inch 


M=s@ bm = 2,0 units/inch 


Pie 5-7 XSCALE 
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C. VARIOUS TYPES OF PENDULUM EXAMPLES 
Example 8: 


Linear approximation to the pendulum, g/l = 1 
No Damping 


mil + X = 0 


PLOT PARAMETERS: 


XSCALE = 2.0 YSCALE = 
XCENTER = 0.0 YCENTER = 
XSIZE = 6.0 YSIZE = 


Size of slopes = 0.2 


Number of slopes per inch = 4 


foae INPUT CARDS USED: 
NELSON, HH. G., PENDULUM 


m+ X = 0 
Z. 0 0.0 6. 0 2.0 0.0 8.0 
ne 4 

me X= 0 
0.0 0. 02 10. 0 
2.0 
1 
0.0 0. 02 10.0 
1.0 
l 
0.0 0. 02 10.0 
4.0 
l 
0. 0 0. 02 10.0 
5.5 
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Example 8: X'' + X= 0 





N . S : * = = cay : | 4 
SO eee 


SSSR RN ppc pce eee 


ee eS eo non Par ao eae ie ee el 0k eee oe 
| 


i 
Cc 


ieee 5-6 mee SCALE — 2.0 units/inch 


iMeoe@ 3B = 2,0 units/inch 
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Example 9: 


Pendulum without damping, g/1 = 1 


X'' + SIN(X) = 0 


PeOT PARAMETERS: 


XSCALE = Ze 
XCENTER = 0. 
XSIZE = 8. 


Size of slopes = 0.2 


Number of slopes per inch = 4 


ima MWNPUT CARDS USED: 


MELSON, HH. G. 
X!'' + SIN(X) = 0 
2, 0 0.0 8. 0 
0.2 4 
om" + SIN(X) = 0 


0 YSCALE 
0 YCENTER 
0 Volz 
PENDULUM 

2.0 0. 0 


0. 0 O02 10.0 
-6. 0 E. 0 

1 

0. 0 0. 02 ro20 
-6.0 2. 0 

] 

0. 0 O02 FO. 0 
-6.0 2.0 

] 

00 OQ. 02 10.0 
-6. 0 4.0 

1 

0.0 0. 02 10.0 
0.0 1.0 

] 

0. 0 we OZ 10. 0 
6.0 1.0 

] 

0. 0 0. 02 LS,.0 
6.0 -2.0 

I 

0. 0 0. 02 10.0 
6.0 -3.0 

] : 
0.0 0. 02 Terre 
6.0 -4.0 
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ee eee eigen ee | gee 


er ae de 


— ae teem Be 


yout/sziun 9°22 = ATVOS-A 
yout/sziun 9°2 = ATIVOS-X 6-G “81g 


te me OE ae wt - -— ae ——_= ie ~ eee Oe ee 


— —, — . ee SO es eae ae ¢ — oe —_ 





0 =(X)NIS + 1X ‘6 oTduex” 
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Example 10: 
Pendulum with velocity damping. 
X'! + A*X'+ SIN(X) = 0 
WHERE: 
A = 0.3 


EPEOLT PARAMETERS: 


BOGALE = Zn YSCALE 
XCENTER = 0.0 MCE INTER 
XSIZE = 8. 0 Work 
Size of slopes = 0.1 
Number of slopes per inch = 4 
met INPUT CARDS USED: 
NELSON, H. G. PENDULUM 
X'' + A*X' + SIN(X) = 0 
Zev) 0. 0 ome Z. 0 0. 0 
0.1 4 
X'' + A*X!' + SIN(X) = 0 
0.3 
0.0 0.05 Zoe0 
-7.0 370 
1 
0. 0 0.05 25,0 
-7.0 5. 0 
1 
0.0 0.05 25.0 
=i. ZnO 
1 
OO 0.05 Zo, 0 
be 0 -5. 0 


70 


“4 It 


a0) 





—— 


yout/sziun 9 °z 


1] 
[x] 
= 
< 
O 
a 
on 


Yout/sztun 9 °z 


tI 
fz] 
= 
<i 
O 
7 
me 
oO 
re 
f 
Ln 
eT 
fx, 








ete -f- Si ae) 


= a — —, 


0 = (X)NIS + .X*V + 4X ‘oT etduexq 
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Example ll: 
Pendulum with velocity squared damping. 
a + ASK *ABS(X') + SIN(X) = 0 
WHERE: 
A=1.0 


EPEOT PARAMETERS: 


ZASCALE = 2.0 YoCALE 
XACENTER = O70 Y CEN SR 
XSIZE = S00 MolZ 


Size of slopes = 0.2 


Number of slopes per inch = 4 


aoe INPUr CARDS USED: 


NELSON, Fi. G. PENDULUM 
me AXA BS(X') + SIN(X) = 0 
2.0 0.0 ono Eaa0 0.0 
0. 2 4 
a + Ae ABS) + SIN(X)-= 0 
1.0 
O:0 G..08 20750 
i. O 4.0 
1 
0..0 O03 20, 0 
-4. 0 4.0 
1 
0.0 O03 20.0 
-7.0 4.0 
1 
0.0 C0 20,0 
a. 0 -4. 0 
] 
0. 0 0. 03 20,0 
ee 0 4.0 
] 
0.0 0. 03 ZO. 0 
a 5 -4.0 


2 





i 


ee eye ee eee —_ — 


Se 7 = ee 
ae ~ 


\ 
\ 
\ 
\ 
\ 
\ 
=e 







a 





X'! + AxXX OKA BS(X!') + SIN( 


X-SCALE = 2.0 units /inch 


ae 


\ 
| 


Example lI: 
Fig. 5-l11 
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Y-SCALE = 220 unitts/imen 





D. SATURATED SYSTEMS 
Example 12: 
A saturated servo system. Saturation occurs when absolute error 


signal exceeds 0.2. 


meets. -0.2) THEN X*' + 0, 2*X' + X = 0.2 
mee Gr, 0.2) THEN X'' + 0.2*xX' +X = -0.2 
me + 0,2*X' + X = 0.0 
PLOT PARAMETERS: 
XSCALE = 0, 3 MS CALE : 0.2 
XCENTER = 0. 0 YCENTER = 0.0 
XSIZE = 6.0 YSIZE = one 
Size of slopes = 0.2 
Number of slopes per inch = 4 
ees INPUT CARDS USED: 
MawSON, H. G. 0902 
SATURATED SERVO SYSTEM 
0.3 0.0 6; 0 OF 2 0. 0 8. 0 
rz + 
ec tal, 0-2, 0)"0,2) THEN X'' + 0.2*X' +X = 0.2 
meee ieee) een x’ + 0.2*x' + X = -0.2 
Pare UneeK + X= 0.0 
0.0 0:05 30.0 
0.7 


74 





Example 12: Saturated Servo System 





2 
e3 
oS 
= — ~~ —_—, ee — — .- — a——ngeei ~<a = 
wee, — ee —— ee oe nee ———— ee 
—. TS, See reg, ee ——_ es 
ea ees — a a ee re 


alleen ys a a 


—_ —— _— Sa eg 
= eg ee ee Ue CO ee 


Fig. 5-12 X-SCALE = 0.3 units/inch 


Y-SCALE = 0.2 units/inch 


(5 





e. 


IDEAL RELAY 
sor. 3x! + 1, 54SIGN(X) = 0.0 


RELAY SYSTEMS 


Example 13: 


PLOT PARAMETERS: 


XSCALE 


XCENTER 


XSIZE 


Size of slopes = 0.2 


Number of slopes per inch = 4 


THE INPUT CARDS USED: 


his 


NMEMSON, H. G. 


0 


©. Zz 


Peet Osow ee + 1. S*SIGN(X) = 0.0 


N © 


Co © 


ie AK RELAY 


0.9 
4. 


6. 


0 


OD: 


0 


0.6 


YSCALE 
YCENTER 


YSIZE 
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= 0.0 


X'! + 0. 3*X' + 1. 5*SIGN(X) 


Example 13: 


i 


; 74 


ee ee ee 
ue 


ae - 
= : 
v ‘ 


o - 


pee, 
PA oe 
a7 oe 


ee ae A 
a f. 





ae ee 
Ze ae ee 


4 





X-SCALE = 1.0 units/inch 


Pay, 503 


Y-SCALE = 0.6 units /inch 


fey 





Example 14: 
Ideal Relay with Rotated Switching Line 
Peet Ov Zoe +), 025I1GN(4.0"X - X*) = 0.0 


meot PARAMETERS: 


XSCALE = fag 0, YOCuEE = 
XCENTER = 0.0 YCENTER = 
2S) VAD = 6.0 YSIZE = 


Size of slopes = 0.2 


Number of slopes per inch = 4 


fae INPUT CARDS USED: 


NELSON, H. G. 

IDEAL RELAY WITH ROTATED SWITCHING LINE 
2.0 0.0 6. 0 a) 0.0 8.0 
0,7 4 

X'' + 0, 2*xX' + 1, OXSIGN(4. O*X - X') = 0.0 


al 80. 0 


See 
NO 
° 
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Beaeople 14: X'' + 0. 2*X' + 1. OFSIGN(4. O¥X - X") = 0.0 






Den a ae SO = Mh 
~ 


aS 
er ra ee r 
ea ee 


a 


ee ee ee ee is ae iS \ 


ae Ne 


a et 
SS SS ee 
ae pe ee ee i i a 
eS a a Pa” ewe ane ae 


Fig. 5-14 X-SCALE = 2. 0 units/inch 


Y-SCALE = 1.0 units/inch 


io 





Example 15: 
Ideal Relay with Rotated Switching Line 
co Oe2 x" + 1 OFSIGN(4. 0*X + X') = 0.0 


PLOT PARAMETERS: 


XSCALE - 2.0 YS CALE = 
XCENTER = OE WMCENTER = 
SIZ. = 6.0 YoiZiz = 


Size of slopes = 0.2 


Number of slopes per inch = 4 


mae INPUT CARDS USED: 


NELSON, H. G. 

Hees RELAY WITH ROTATED SWITCHINGHINE 
2.0 0.0 6. 0 Cs: Oe 8.0 
0.2 a 

mea 8, 2x + 1 Os SIGN(4. OFX +X") = 0.6 


ol 20. 0 


U1 O © 
U1 O © 
© 
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Peete x + 1. OFSIGN(4. OFX + XX") = 0.0 


Example 15: 


a 


oe ae ee 
= ee a ae ae 


aa ae ae 
eae ae eae 
ae eae 


as Sepa | 


ae oo eee eee 
Pom ol Mee 
ae eS ae 





go ee ee” ae eee 


; 4 ye a 


x Te ae 


ee ge ee a ee ee ae 





Vee eee 


Reamer = 2. 0 units/inch 


ire. 5=15 


Y-SCALE = 0. 8 units/inch 
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Example 16: 
Ideal Relay with Rotated Switching Line 
pero Over + 1, O25IGN(2Z.0sX + X') = 0.0 


meor PARAMETERS: 


XSCALE = Z2.0 YSCALE = 
XCENTER = 0.0 YCENTER- «= 
XSIZE = 6. © hora = 


Size of slopes = 0.2 


Number of slopes per inch = 4 


aa: INPUT CARDS USED: 


NELSON, H. G. 

Peel RETA Y WIitH ROTATED SWITCHING, LINE 
2.0 O20 G0 i 0 Ome &..0 
0.2 4 

Bo OL 25x! + 1, 045IGN(Z2. OFX + X') = 0.0 


20. 0 


mm Oo © 
Mm © © 
© 
—_ 
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meotaple to; x ' + 0. 2*x' + 1. OFSIGN(Z. O¥X + X") = 0.0 









ci hecceecetan~ tests 
t VE a 
ea) ee, ee 


CL oe ew —_— ~~ . S 





= pe pr gerne. oe = ss clam 


OTDM Me ee Se 
Vy ae a 
i 7 ee ae 


ee OS. \ 


oa cas; 
t2 


SN ™N 


Sie 7 oy aes 
Ue ar A Ye 


| -_ -_ — 


On yet ee 


— 


= Ss 
Peek SS 
\ 
\ 


mode ee 


~~ 


Powe to Meo CALE, = 2. 0 units /inch 


Y-SCALE = 1.0 units /inch 
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Example 17: 
Ideal. Relay with Rotated Switching Line 
pet One «xX + 1, OF5IGN(Z. OFX - X') = 0.0 


Peeler aRAMETERS: 


XSCALE - ZO YoerlLE 
SemNIER = 0.0 WGENT ER 
XSIZE = 6.0 YSIZE 


Size of slopes = 0.2 


Number of slopes per inch = 4 


PiteNeOor CARDS USED: 


MELSON, HH. Cc. 


IOEAL RELAY WITH ROTATED SWITCHING LINE 


e. 0 0.0 6.0 0.8 0. 0 
Ore = 
eee One ee + 1 OSIGN(2Z, OFX - X') = 0.0 


oul 80. 0 


moo 
mn © © 
Oo 
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cee One x + 1. OFSIGN(Z. OFX =~ X") = 0.0 


Example 17: 


a ee 
| 
ae = weal 


ea ee ee ee 


eee ee 


Se or ee eee aie 
a 


eel 





s 
Sora ad 


ao oe 


Fig. 5-17 X-SCALE = 2.0 units/inch 


Y-SCALE = 0. 8 units/inch 
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FF. DEAD ZONE 
Example 18: 


iad, -A,OR,X.GT,A) THEN X'' + B¥X' + C¥X = CXA4SIGN(X) 
X'' + B¥X' = 0.0 


WoeR E: 
A= 0.3 
B= 0.2 
ee 1.0 


PLOT PARAMETERS: 


XSCALE = 0.5 YSCALE = 0.3 
XCENTER = O20 YCENTER = 0.0 
mSIiZ = 620 YSIZE = 8.0 
Size of slopes = 0.2 
Number of slopes per inch = 4 
ima INPUT CARDS USED: ° 

NELSON, H. G. 

DEAD ZONE 
0.5 0.8 6.0 0.3 0.0 Oued 


ee 4 
ae Lt. -A. ORE XGA) THEN X" + BX! + C¥X = C™A*SIGN(X) 
X'' + BeX' = 0.0 


G23 OFZ 1.0 
OO 0.1 507 0 
Gr. © i 

] 

O20 0.1 50. 0 
EO 
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Example 18: Dead Zone 

Myeews. 2, Ob] xX, GT, A) THEN X"' + B¥X' + C*#X = GC¥A*SIGN(X) 
peat B=X'e= O00 

m= 0. oe be Oud, C = 1.0 


/ es ei a a x x \ x 
ee NS \ \ 

; s ~ = A SS \ . 
/ A 
/ a 
u \ \ \ 


~ 
— 


- - 
© 


——— — 
a 
— 


-" 
— - 


— 


Fae a ph ptf 


— 
-_ 


—— ee 


---+--+h ry Lb, 





Fig. 5-18 X-SCALE = 0.5 units/inch 


vene@ oe = 0,3 units /inch 
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Example 19: Dead Zone 


Meeso( So), LT. 2.0) THEN X"' + 0.2*X' = 0.0 
pe + .05 2x" +t X = 0.0 


PLOT PARAMETERS: 


XSCALE ~ a0 YSCALE - 
XCENTER = 0.0 YCENTER - 
XSIZE = 8. 0 Mow = 
Size of slopes = 0.2 
Number of slopes per inch = 4 
THE INPUT CARDS USED: 

NELSON, H. G. 
DEAD ZONE 
z@. 0 0.0 8.0 Zao 0. 0 Oan0 


w.2 4 
Trea bes). LT. 2. 0)sTHEN X'' + 0. 2*X' = 0.0 
pe 7+ OR 2=xX' + X = 0.0 


O05 a). 0 


O& © © 
OO © 
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IF(ABS(X). LT. 2.0) SHEN KX" + 0lZee ee CRO 


X' +0, 248? +X = 


Example 19: 


0. 0 


aie SO A a ae 
| 


f 





i aS Se 
a 1 Oe 
J 


N 


7 oS) 
| | 


4 ~~ 
~~ eo. = a a 
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ou 


Fig. 5-19 X=SCALE = 2.0 umitey inch 


Y-SCALE = 2. 0 units/inch 





G. A SECOND ORDER LINEAR SYSTEM 
Example 20: 
The following three plots demonstrating the effect of various 
damping were the result of the input cards listed below. 
m+ 2. O*A*B eX! + B*¥2*X = 0. 0 
Where: 


A 
B 


is the natural frequency 
is the damping factor 


PLOT PARAMETERS: 


XSCALE - 2.0 YSCALE = fap’) 
XCENTER = 020 YCENTER = Oe 
XSIZE = 8. 0 YSIZE = 6.0 


Size of slopes = 0.2 


Number of slopes per inch = 4 


mab INPUT CARDS USED: 


NELSON, H. G. 
X'' + 2. O¥A*¥B¥*X' + B¥*2*X = 0.0 


ZU 0.0 Gal) Zr) O70 0.0 
OE 2 = 
bet 2, 07A+BsX' + 6¥%*2*X = 0.0 
1.0 130 
0. 0 rae 10. 0 
( O 
] 
OO Oe £O.0 
-7. 0 
] 
O.0 0.1 20 
3205 
IL 
0. 0 0.1 BO? 0 
7s. 
2 
G5 1,0 
O90 0.1 15.0 
1.0 
] 
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0.0 


X!'' + 2, OXA*B*X! + Bowe EX 


A=1.0, B=1I1.0 


Example 20(a) 


ee 


ie) | 


NS neo 


i 


ee gn 


a pe 
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a i 





2.0 units/inch 


XA-SCALE 


Fig. 5-20(a) 


Y-SCALE = 2.0 units/inch 





0. 0 


X'! + 2, OFA*BXX! + B*KZ XK 


Example 20(b) 















SD 
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e a 
ae ee 
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vo 


mah Uc 
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Zz A Caz 
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Pe r 
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ee ee 


— 
Cc? 
fame 
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See eee oo cee 
Re ee ee 
eet OO 


<2 
| 
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es) ee eee a> eet ee 


ee OO a 


-S 


X-SCALE = 2.0 units/inch 


Fig. 5-20(b) 


2.0 units/inch 


Y -SCALE 





0.0 


X!" + 2, OXA*B*X! + BxKZ xX 


Example 20(c) 


ae 
vo 
a 


Se ns 
\ » a ee ee 
eS. ~~. a ~~ HH eae 
ee ee Seepere ys 
x | X = _ = ee nt ae 
oe NNR ROS eee peor ae ao a ie 
—— a ~~ >_> — —— LO Le yee 
ad a, eee ay a 


| 
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ee a —_— oe — ——- — o_O Sir oe eR 
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X-SCALE = 2. 0 units /inch 


Fig. 5 -20(c) 


Y-SCALE = 2.0 units /inch 





Pe SeePiIPTH OR DER*S YSTEM 
Example 21: 
mee ST. 0.10015) SKIP 3 
= 3000511.6 +A 
haz = 265120.0 + A*351T. 0 
V3 = A*38500. 0 
eee eee ees S55, SFX"! + VI¥X'' + VZeX' + V3=X = 0.0 
WHERE: 
G = 30212. 0 


PLOT PARAMETERS: 


RoOCALE = Zeno) SGA) = 2020 
XCENTER = Ono YCENTER = 0.0 
mo ZB = 6.0 YoIZ | = 8. 0 


fee INPUT CARDS USED: 


MELSON, H. G, 
mere ORDER SYSTEM . 
Z. O 0.) Gnd ZOOM OO 8. 0 


IF(T.GT. 0.0015) SKIP 3 

Vl = 3000511.6+A 

V2 = 285120. 0 +A¥*3511.0 

V3 = A*38500. 0 

Kt + 122, 2X + 50835, SeK'!' + VEX! + V2%X' + V3X = 0.0 


Boz12. 0 

FO 0.001 0.4 
=n 0 

] 

of 0 U2 00L 0.4 
ZO 

1 

0.0 0.001 0.4 
=5, 1 

] 

O. 0 0.001 0.4 
4.6 


5 





Example 21: A FIFTH ORDER SYSTEM 

IF(T .GT. 0.0015) SKIP 3 

V1 = 3000511.6 +A 

V2 285120. 0 + A*3511,. 0 

V3 = A*38500. 0 

X™ 1 4 122, 2eK'' + 50835, BX"! + VIEX'! + V24K! + V3¥X = 0.0 


cc 
a 
a a a ee dE 





2.0 units /inch 


es 5-21 X=SCGALE 


Y-SCALE = 20. 0 units/inch 
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/ 


1. SHYSPERESE AND DEAD ZONE 

Example 22: 
RELAY Vite DEAD ZONE AND HYSTERESIS 
THE SYSTEM DESCRIPTION EQUATIONS AND/OR PROGRAM: 
V3 =T 
IF(ABS(X) .LE. A)VI = 0.0 
IF(ABS(X) .GT. B)V1 = SIGN(X) 


par + 0,2*X' + V1 = 0.0 


WHERE: 
m= 0. 8 
Beal. 2 


The first statement is used to suppress the slopes. 
miOr PARAMETERS: 


XSCALE = Dae YSCALE = 07> 
XCENTER = 020 YCENTER = 0.0 
ASIZE = 6.0 YSIZE = San 


moe etNPOLT CARDS USED: 


NEESON, H. G. 
meen yY Well DEAD ZONE AND HYSTERISES 


1 0 0. 0 6. 0 0.5 0. 0 8. 0 
V3=T 
[Pees x) bm. A) Vi = 0.0 
IF(ABS(X), GT, B) V1 = SIGN(X) 
X'' + 0,2*X' + V1 = 0.0 
0.8 Ne 
0.0 0. 04 30.0 
22.8 


ot 





Example 22: RELAY WITH DEAD ZONE AND HYSTERESIS 


f 





Ie > -22°° X-SCALE = 1. 0 units/inch 


Y-SCALE ~- 0.5 units/inch 
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J. 


esGOND ORDER SYSTEM WITHA RESTORING FORCE CORRESPONDING 
TO A NON-LINEAR SPRING 


Example 23: 


NON-LINEAR RESTORING FORCES 


Dae SYSTEM DESCRIPTION EQUATION: 


X'' + 0,9*X' + (X + A# 
WHERE: 
02 
PLOT PARAMETERS: 


XSCALE 


ACENTER 


ASIZE 


=— 
— 


I 


Size of slopes = 0.2 


Number of slopes per inch = 4 
foe (INPUT CARDS USED: 
NELSON, H. G. 


SECOND ORDER SYSTEM WITH NON-LINEAR SPRING 
6.0 


ZaaO 
ec 


0.0 


+4 


#3) = 0. 0 


YSCALE 
YCENTER 


YSIZE 


0. 0 


he 0, 92x"! + (Xt AXK%%3) = 0.0 


Oo © He ° 
-o° Bon 
© = © 


wm 
ce 
Oa) 


© 


Mort OorF-torrortioirtwnNioo 
es #6 in °* e e« 6 e & 
5 O © 
N 


CO © 


0402 


0. 02 
-7.0 


G02 
-5. 0 


LS. 


Lo. 


ree 


ioe 


Lon 


hen 


0 


oo 





X'! +0, 9%X! + (X + A*X**3) = 0.0 


Example 23(a): 


0.1 (gives a ''soft'' spring) 


A 


ed 
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S\N Eee ae 
a Oa SN NE SS 
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<3 Sh NS 
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MeSGALE = 2.0 units/inch 


Fig. 5-23(a) 


Y-SCALE = 2.0 units/inch 
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0.2 (gives a "hard'' spring) 


poe Toe + (X + A*X*43) = 0.0 


A 
a ee eo “ 


3 
oe 


Example 23(b): 


ff 
ee ee = 
/ 
/ 
f 


/ i / 
Bo? 
yo fo 
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, ei 





f 
/ 
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MeasG@ wee = 2,0 units/inch 
Vec@etr = 2, 0 units/inch 


Fig. 5-23(b) 





Example 24: 
A FORCED SYSTEM WITH A NON-LINEAR RESTORING FORCE 
MmdteorolEM DESCRIPTION: 


eee ge + (XS Ol 2x3) = A*COS(B*T) 


Maik i: 
me = 3. 0 
pe] 0. 7 


EPeOT PARAMETERS: 


oO Li = 0 YSCALE = tO 
XCENTER = GEO YELNPTER = wo 
XSIZE = or0 WoIZE = 3. 0 


THE INPUT CARDS USED: 

DMMBSON, H.. G., 

FORCED NON-LINEAR SECOND ORDER SYSTEM 
1.0 0. 0 6.0 0 0.0 Sead) 


me St. Jax + (XxX + 0. 2X3) = A*COS(B*T) 


30. 0 


WOOW 
Oo aS © 
Oreo 
Oo ~) 

at 


NEESON, oH. G. 
FORCED NON-LINEAR SECOND ORDER SYSTEM 
ZO Oa ond ZO GeO 8.0 


X'' + (X + 0. 2*¥X*%3) = A*COS(B*T) 


5.0 0. 8 

0.0 0. 04 5 050 
0...0 Ore 

2 

50 0 

0. 0 0. 04 510). © 
0.0 eee 


BZ 





Example 24(a): X'' + 0.9*X' + (X + 0.2*X**3) = A*COS(B*T) 
pees. 0, ©B = 0.7 





Me. Saedia) X-SCALE = 1.0 units/inch 


1.0 units/inch 


Y-SCALE 
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Example 24(b): X'! + 0. O*X' + (X + 0.2*X*%3) = A*COS(B*T) 
A=5.0, B=0.8 


un 
CI 
=) 


-UC6 





awd 


Fig. 5-24(b) X-SCALE = 1.0 units/inch 


Y-SCALE = 1.0 units/inch 
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Example 24(c): X!'' + 0. O*X' + (X + 0.2*X**3) = A*COS(B*T) 
me eo, 6B = 1.0 


Jos 





my 
ct 
a 


Fig. 5-24(c) X-SCALE 


i 


1. 0 units /inch 


Y-SCALE = 1.0 units/inch 
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K. BANG-BANG SYSTEMS 
Example 25. 
A SECOND ORDER BANG-BANG SYSTEM 
erolTEM EQUATION: 
ee!) + SIGINGEX +00n5:*X" *HBS(XK")) =-02-0 


PLOT PARAMETERS: 


XSCALE = 2.0 YSCALE 
XCENTER = 0.0 Yeu NTER 
ASIZE = 8.0 YSIZE 


Size of slopes = 0.2 


Number of slopes per inch = 4 


mmr INPUT CARDS USED: 
BEESON, HH. G. 
SECOND ORDER BANG-BANG 
20 UE U 8. 0 Lo 0. 0 


0.2 x 
pee tT olen X + 0. SX A BS(X")) = 0.0 


70 OMOZ enn 


0 

] 

0.0 0.02 E5220 
120 
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Example 26. 
A THIRD ORDER BANG-BANG SYSTEM 
BYOoTEM EQUATIONS: 
V1 = SIGN(X! + 0. 5*X!'*ABS(X"')) 


ez 


X + (1. /3. RXR + VLEX RX! 
V3 = V1 *(0. SX! 2 + VIX) RL, S 
po + SIGN(VZ + V3) = 0.0 


PLOT PARAMETERS: 


XSCALE = 2.0 YSCALE 
AXACENTER = OO YCENTER 
XSIZE = 6.0 YSIZE 


The input cards used: 
DBSLSON, H. G. 
THIRD ORDER BANG-BANG SYSTEM 


a, 0 0. 0 6.0 ioe 0. 0 


lenge 1.0, 54x =A BS(X'')) 
Pere a Se ees + VIX! tae 
V1*(0. 5B&eX ke? + VLkX'! yaa] , 5 
me + SIGN(VZ + V3) = 0.0 


< 
IN 
Wow 
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A THIRD ORDER BANG-BANG SYSTEM 


Example 26: 





GZS 


Fig. 5-26 X-SCALE = 2. 0 units/inch 


1. 0 units /inch 


Y-s@2bE. 
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Pee Vemioreons a ND SUGGESTIONS FOR FURTHER RESEARCH 


This thesis has introduced an easy to use, flexible, and powerful 
analysis tool for the study of non-linear ordinary differential equations. 
The primary application emphasis has been in relation to control 
system problems for which many examples have been presented. It 
has been shown that the use of an automatic built-in compiler coupled 
with interactive graphics result ina very useful analysis program. 

For a large number of quite diverse examples it has been demonstrated 
in this thesis that one can very adequately specify a problem using 
data cards or by typing on an interactive graphics screen. 

This study has suggested several basic ideas for further research. 

1. An extension of the program following the development presented 
in Chapter II would give the program the capability of solving systems 
described by several simultaneous non-linear differential equations. 

2. If one were concerned with execution efficiency, it would be 
very worthwhile to rewrite the interpreter in assembly language. The 
interpreter is called four times for each solution point and once for 
each slope line. Thus the slope generation and solution generation 
times could be speeded up considerably with a faster interpreter. The 
other parts of the program are used relatively infrequently and there- 
fore efficiency is not an important factor. 

3. It would be worthwhile to provide a choice of integration 
schemes to the user in addition to the present fourth order Runge- 


Kutta. 





4. A built-in compiler could be used in several other applications. 
One of the very obvious areas is that of optimization. For example, 
avery useful directed search algorithm called ''DIRECT"' is available 
on the XDS 9300. The function to be minimized is incorporated into a 
function subroutine which must be supplied by the ee A modification 
of the compiler and interpreter presented in this thesis could be used 
as avery powerful interface for describing the function and its 
boundaries. 

5. It would be very worthwhile developing methods to assist in 
the analysis of systems higher than second order. In general, the 
presentation of the grid of slopes is not worthwhile for systems of 
higher than second order since slope lines are ina space which has a 


dimension equal to the order of the differential equation. 





- APPENDIX A 


PeOUIRED JOB CONTROL LANGUAGE FOR BATCH VERSION 


1. When using the compiled binary object deck use the following 
procedure: 


Vi) JOB CARD -~Green- 
// EXEC FORTLG, REGION. GO=160K 
//LINK.SYSIN DD * 

BINARY OBJECT DECK 
//GO,.FTO6F001 DD SYSOUT=A,SPACE=(CYL, (6)) 
//GO.SYSPLOT DD SYSOUT=C, SPACE=(TRK, (1, 10)) 


//GO.SYSPLOTS DD UNIT=S YSDA, SPACE=(TRK, (2, 6)) 
//GO.SYSIN DD * 


DATA 
[* ~Orange- 
2. If using the source deck use the following procedure: 

// JOB CARD -Green- 
// EXEC FORTCLGP, REGION. FORT=150K, REGION. GO=160K 
//FORT.SYSPRINT DD SYSOUT=A,SPACE=(CYL, (1, 1)) 
//FORT.SYSIN DD * 

FORTRAN SOURCE PROGRAM 
/% 
//GO. FTO6F001 DD SYSOUT=A, SPACE(CYL, (6)) 
//GO.SYSPLOT DD SYSOUT=C, SPACE=(TRK, (1, 10)) 
//GO.SYSPLOTS DD UNIT=SYSDA,SPACE=(TRK, (2, 6)) 
//GO.SYSIN DD * 

DATA 


/* Orie e = 


Pr 


— 





i APPENDIX B 


1. Required control cards for Interactive Graphics Version. 


BOOT CARD 


Pal CH DEGK 


46 AGT 
46 JOB 
& ASSIGN BI =MTOA 
4 ASSIGN GO =MTOA 


& LOAD XM, MAP 
O&O SEG MAIN-(SINIT, SSETUP2, SCOMP, (SGRID-SINTER), (SSOLVE- 
& DATA SINTER )) 
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PReoocRkAM A 


at 


LISTING OF BATCH VERSION 
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